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Epigenetic changes and functional study of 
HOXAT11 in human gastric cancer 


Aim: To examine epigenetic changes and the function of HOXA11 in human gastric 
cancer (GC). Materials & methods: Seven GC cell lines, five cases of normal gastric 
mucosa and 112 cases primary GC samples were used in this study. Results: Expression 
of HOXAT1 and lack of promoter region methylation were found in NCI-N87, MKN45, 
BGC823 and HGC27 cells. Loss of expression and complete methylation were found 
in AGS gastric cancer cells. Reduced expression and partial methylation were found 
in MGC803 and SGC7901 cells. Restoration of HOXA11 expression was induced by 
5-aza-2'-deoxycytidine. HOXA11 was methylated in 81.25% (91/112) of primary GCs. 
The presence of methylation was associated with male gender, tumor size, tumor 
differentiation and lymph node metastasis (all p < 0.05). Restoration of HOXAT1 
expression reduced cell proliferation, invasion, migration and induced apoptosis and 
G2/M phase arrest. HOXA11 was found to inhibit Wnt signaling by upregulating NKD1 
expression. Conclusion: Epigenetic silencing of HOXA11 promotes GC proliferation, 
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Gastric cancer (GC) is the fourth most 
common malignancy and ranks the second 
leading cause of cancer-related death in 
the world [1,2]. Salt, nitrite, alcohol, smok- 
ing and Helicobacter pylori infection were 
regarded as important risk factors [3]. Even 
though aberrant genetic changes play criti- 
cal roles, they are not enough to explain the 
mechanism of gastric carcinogenesis [4,5]. 
Epigenetics is becoming one of the impor- 
tant fields of cancer research. Increasing evi- 
dence suggests that epigenetics is involved 
in gastric carcinogenesis [6,7]. HOX genes 
are evolutionarily conserved genes encoding 
transcription factors that regulate mamma- 
lian embryonic development [8,9]. HOXAII 
plays an important role in regulating cell 
differentiation and proliferation [10,11]. Pro- 
moter region hypermethylation of HOXAI1 
was found in different human cancers [12-15]. 
The epigenetic changes and the function of 
HOXAII in human gastric cancer remain 


unclear. In this study, we analyzed the epi- 
genetic regulation and the mechanism of 
HOXAII in the development of human 


gastric cancer. 


Materials & methods 

Primary GC samples 

One hundred and twelve cases of GC sam- 
ples were collected from the Chinese PLA 
General Hospital. The median age of can- 
cer patients was 58.5 years (17—85 years). 
Tumor stage was classified according to 
TNM staging (NCCN2012), including 
stage IIA (n = 4), stage IIB (n = 12), stage 
IITA (n = 9), stage IIIB (n = 48), stage IIIC 
(n = 15), stage IV (n = 24). Five cases of 
normal gastric mucosa were collected from 
noncancerous patients. All samples were 
collected according to the approved guide- 
lines of the Chinese PLA General Hospi- 
tal’s institutional review board and kept 


at -80°C. 
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Cell culture & 5-aza-2'-deoxycytidine treatment 
Seven GC cell lines were included in this study (AGS, 
NCI-N87, SGC7901, MGC803, BGC823, HGC27 
and MKN45). The GC cell lines were established from 
primary human gastric adenocarcinoma cells previ- 
ously and maintained in 90% RPMI 1640 (Invitro- 
gen, CA, USA), supplemented with 10% fetal bovine 
serum and antibiotics [16]. GC cells were split to a low 
density (30% confluence) 12 h before treatment, then 
treated with 5-aza-2'-deoxycytidine (5-aza; Sigma, 
MO, USA) at a concentration of 2 uM in the growth 
medium, which was exchanged every 24 h for total 
96 h treatment. At the end of the treatment course, 
RNA was isolated as described below. 


RNA isolation & semiquantitative RT-PCR 

Total RNA was isolated by Trizol reagent (Invitro- 
gen, Carlsbad, USA). RNA quality and quantity were 
evaluated by 1% Agarose gel electrophoresis and spec- 
trophotometric analysis. Semiquantitative reverse tran- 
scription-PCR (RT-PCR) was performed as described 
previously [17]. RT-PCR primers were as follows: (F) 
5'-GCAGCAGAGGAGAAAGAGCG-3' and (R) 
5'- GTGGATT TGCTGAGTAGTACTG-3". 


DNA extraction, methylation-specific PCR & 
bisulfite sequencing 

Genomic DNA from GC cell lines and GC tissue 
samples was prepared by proteinase-K method. Meth- 
ylation-specific PCR (MSP) and bisulfite sequenc- 
ing (BSSQ) were performed as described previously 
[18,19]. To obtain optimal conversion of cytosine to 
uracil by bisulfite treating, 1 ug genomic DNA was 
used, because increasing the amount of input DNA 
carried the risk of incomplete conversion and led to 
false-positive results [20]. Sodium bisulfite treated- 
DNA should be treated as RNA and stored at-20°C to 
prevent further degradation. It is important to avoid 
repeated freezing and thawing. Bisulfite sequencing 
was used to validate MSP result and provide complete 
CpG methylation information on all CpG-sites within 
the PCR product. 

Each MSP reaction incorporated approximately 100 
ng of bisulfite-treated DNA, 25 pmoles of each primer, 
100 pmoles dNTPs, 2.5 ul 10x PCR buffer and 1 unit 
of JumpStart Red Taq Polymerase (Sigma) in a final 
reaction volume of 25 pl. Cycle conditions were: 95° 
C x 5 min; 35 cycles x (95°C x 30 s, 60°C x 30 s, 
72°C x 30 s); 72° C x 5 min. MSP products were 
analyzed using 2% agarose gel electrophoresis. To 
interpret the MSP results properly, positive, negative 
and H2O controls were included. MSP primers and 
BSSQ primers were designed around the transcrip- 
tional start site. To accomplish optimal discrimina- 


tive power between methylated and unmethylated 
DNA, the primers should have a substantial num- 
ber of non-CpG cytosines in the original template 
and contain one to three CpG dinucleotides in the 
3 region of the primers [20]. The sequences of HOXAII 
MSP primers were as follows: (UF) 5’TGTGTGT- 
GAAGTGATTTTTAGAGAGTAT-3’, (UR)5’-AAC 
AATCTCTATACACAAACTCCTCCA-3'; 
(MF)5'-TGCGCGAAGTGATTT 
TTAGAGAGTAC-3’, (MR)5’-CAATCTCTATACAC 
GAACTCCTCCG-3". Primersequences of BSSQ were: 
(F)5'-GTTTYGGGTTTAGATTTTTTTAGT-3', 
(R)5'-TAACCRAACTCTTAACCAACAC-3’. The 
unmethylated (U) primer will only amplify sodium 
bisulfite converted DNA in unmethylated condi- 
tion, while the methylated (M) primer is specific for 
sodium bisulfite converted methylated DNA. Sample 
was regarded as unmethylation when PCR product 
was amplified only in unmethylation reaction and 
sample was regarded as methylation when amplified 
only in methylation reaction. Sample amplified in 
both unmethylation and methylation reaction was 
regarded as partial methylation. Since clinical tissue 
samples should be considered heterogeneous, both 
unmethylation and methylation amplicons can be 
expected. 


Immunohistochemistry 

Immunohistochemistry (IHC) staining was per- 
formed in 45 cases of available matched cancer and 
adjacent noncancerous tissue samples. The procedure 
was performed as described previously [17]. Immuno- 
histochemistry (IHC) was performed on 4 um thick 
serial sections derived from formaldehyde-fixed paraf- 
fin blocks of colorectal cancer and paired adjacent tis- 
sue. After deparaffinization and rehydration, endog- 
enous peroxidase activity was blocked for 30 min in 
methanol containing 0.3% hydrogen peroxide. After 
antigen retrieval, a cooling-off period of 20 min pre- 
ceded the incubation with the primary antibody. 
Anti-HOXA11 Mouse Monoclonal (8B8) was diluted 
to 1/200 (LifeSpan BioScience, MD, USA), anti-p- 
B-catenin was diluted to 1/100 (Bioword, Shanghai, 
China). The staining intensity and extent of the 
staining area were graded according to the German 
semiquantitative scoring system as described before 
[1721,22]. Staining intensity of the nucleus, cytoplasm 
and/or membrane was characterized as follows: no 
staining = 0, weak staining = 1, moderate staining = 2, 
strong staining = 3; the extent of staining was defined 
as: 0% = 0, 1-24% = 1, 25-49% = 2, 50-74% = 3, 
75-100% = 4. The final immune-reactive score 
(0-12) was determined by multiplying intensity score 
by the extent of staining score. 


202 


Epigenomics (2015) 7(2) 


fsg 


future science group 


HOXA11 expression plasmid construction 

The full length complementary DNA (cDNA) of 
HOXAII was obtained by RT-PCR. The primers were as 
follows: HOXAII-EGFPNI-F: 5'- C@GAATTCATG- 
GATTTTGATGAGCGTGGTCCCT-3'; HOXALI- 
EGFPNI-R: 5'- CGGGATCCCGGAGGAGTG- 
GATTTGCTGAGTAG-3’.. The products were 


subcloned into peDNA™ 3.1 expression vector. 


shRNA knock down assay 

ShRNA sequences of HOXAII are as follows: (F) 5'- 
GATCCGCCCAATGACATACTCCTACTTTC 
AAGAGAAGTAGGAGTATGT- 
CATTGGGCTTTTTTGGAAA-3'; (R) 5'-AGCTT 
TTCCAAAAAA GCCCAATGACATACTCCTAC- 
TAAGTTCTCTAGTAGGAGTATGT- 
CATTGGGCG-3’.. The scramble control non- 
targeting short-hairpin RNA sequence was 
5'-TTCTCCGAACGTGTCACGT-3’ (Ambion Inc., 
TX, USA). 


Transfection 

AGS cells were transfected with pCDNA3.1-HOXAI1 
or the empty vector using FuGENE 6 (Roche Applied 
Science, IN, USA) according to manufacturer’s instruc- 
tions. BGC823 cells were transfected with HOXAII- 
shRNA or scramble control using FuGENE 6 (Roche 
Applied Science). 


Colony formation & cell viability detection 

Cells were seeded at a density of 2000 cells/well in 
6-well plates to grow for 2 weeks. Then cells were fixed 
with 75% ethanol for 30 min, stained with 0.2% crys- 
tal violet (Beyotime, Nanjing, China) for 20 min and 
then the number of clones was counted. Cells were 
plated in 96-well plates at a density of 2 x 10° cells/well, 
and cell viability was measured at 24, 48 and 72 h 
using the MTT assay (KeyGEN Biotech, Nanjing, 
China) according to the manufacturer’s instruction. 
Absorbance was measured on a microplate reader 
(Thermo Multiskan MK3, MA, USA) at a wavelength 
of 490 nm. 


Transwell assay for migration 

The effect of HOXAII on cell migration was detected 
by using Transwell Assay System (Corning, High 
Wyombe, UK). Cells were harvested and suspended 
in the serum-free medium and then placed into the 
upper well at a concentration of 2 x 10° cells/200 ul. 
The complete medium with 20% fetal bovine serum 
was placed into the lower well (600 ul). The chamber 
was incubated for 24 h. The cells still on the upper 
surface were scraped gently and washed out with 
PBS for three-times. The cells migrated to the lower 
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surface of the membrane were stained with crystal 
violet and counted in three independent high-power 


fields (x200). 


Transwell assay for invasion 

A total of 1 x 10° cells were suspended in 200 ul of 
serum-free medium and loaded onto the upper compart- 
ment of an invasion chamber containing a polycarbon- 
ate membrane with an 8 um pore size which was coated 
with a layer of extracellular matrix (ECM; Matrigel™, 
BD, NJ, USA). After 24 h of incubation, the invasive 
cells migrating through the ECM layer to the com- 
plete medium in the lower compartment were stained 
with crystal violet and the number of invasive cells was 
counted in three independent high-power fields (x200). 


Flow cytometry analysis 

Cell phase distribution was examined as described 
previously [23]. For apoptosis analysis, Annexin 
V-FITC/PI Apoptosis Detection Kit (KeyGEN Bio- 
tech) was conducted according to manufacturer’s 
instructions. 


Gene expression array analysis 
RNA from HOXAII unexpressed and re-expressed 
AGS cells was purified using the RNeasy Mini Kit 
(QIAGEN, Shanghai, China) and the expression was 
analyzed by expression array (Agilent Human 4 x 44k, 
Bio Corporation, Beijing, China). 


Dual-Luciferase reporter assay 

AGS cells were seeded at 2 x 10*/well in 96-well 
plates and incubated for 24 h. The method was 
described previously [17]. To examine transcrip- 
tional activity driven by TCF/LEF, AGS cells were 
transfected with 20 ng/well Topflash reporter vector 
(TCF/LEF-responsive reporter, containing 4x TCF- 
binding sites); 2 ng/well pRL-TK and pCl-neo-B- 
catenin expressing vector (20 ng/well) were used to 
activate the reporter gene. Asa negative control, another 
group of AGS cells were transfected with 20 ng/well 
Fopflash reporter vector (containing mutant TCF/LEF 
binding sites), 2 ng/well pRL-TK control vector and 
B-catenin expressing vector (20 ng/well). To detect 
the effect of HOXAII on B-catenin/TCF luciferase 
reporter activity, pEGFP-HOXAJ/ and wild-type or 
mutant B-catenin were cotransfected into AGS cell. 
Forty-eight hours after transfection, relative lucifer- 
ase activities were measured with the Dual Luciferase 
Reporter Assay System (Promega, Shanghai, China) 
according to the manufacturer’s protocol. Each experi- 
ment was repeated for three-times. Statistical analy- 
sis was carried out by Student’s t-test, p < 0.05 was 
considered as statistically significant difference. 
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Figure 1. HOXA11 expression was regulated by promoter region methylation in gastric cancer (see facing page). 
(A) The expression of HOXA11 was detected by semiquantitative RT-PCR before and after 5-aza-2'-deoxycytidine 
(5-AZA) treatment in gastric cancer cell lines. GAPDH: internal control. (+): 5-AZA treatment, (-): no 5-AZA 
treatment; AGS, NCI-N87, SGC7901, MGC803, BGC823, HGC27 and MKN45 are gastric cancer cell lines. (B) Upper: 
HOXA11 methylation status in gastric cancer cell lines. Lower: Bisulfite sequencing of HOXA11 promoter region 
(+409 bp to +135 bp) in AGS, BGC823 and SGC7901 cell lines. MSP: MSP amplification region, spanning 274 bp. 
Filled circle: methylated CpG site; open circle: unmethylated CpG site. (C) Upper: Representative MSP results of 
HOXAT1 in human gastric cancer. Lower: HOXA11 MSP results in. normal gastric mucosa. (D) Left: IHC results of 
HOXAT1 expression in gastric cancer and adjacent tissue samples (200x, upper), and the expression of HOXA11 
in promoter region unmethylated and hypermethylated gastric cancer samples (200x, lower). Right: RT-PCR and 
western blot results in gastric cancer and normal gastric mucosa. 

C: Gastric cancer; IVD: /n vitro methylated DNA served as methylation control; M: Methylated alleles; N: Normal 
gastric mucosa; NL: Normal peripheral lymphocytes DNA served as unmethylated control; U: Unmethylated alleles; 


H,O: Water. 


Western blot 

Western blotting was performed as described previ- 
ously [17]. Antibodies were diluted according to manu- 
facturer’s instruction. Antibodies employed were as 
follows: antic HOXA11 Mouse Monoclonal (8B8) Anti- 
body (LifeSpan BioScience, MD, USA), rabbit anti- 
Cylin D1 (Bioworld Tech., MN, USA), rabbit anti-c- 
Myc (Bioworld Tech.), rabbit anti-Survivin (Bioworld 
Tech.), rabbit anti-mmp2 (Bioworld Tech.), rabbit anti- 
mmp9 (Bioworld Tech.), rabbit anti-cyclinB1 (Bioworld 
Tech.), rabbit anti-B-catenin (Epitomics, CA, USA) and 
rabbit anti-p-B-catenin (CST, MA, USA). 


Statistical analysis 

SPSS 17.0 software was employed. y? test and Student’s 
t-test were used. Quantitative data were presented 
as the mean + SD. A p-value of less than 0.05 was 
considered as statistical significance. 


Results 

HOXAT1 expression is controlled by promoter 
region methylation in human GC cell lines 

To examine the expression of HOXALI, semiquantita- 
tive RT-PCR was employed. Expression of HOXAII 
was found in NCI-N87, MKN45, BGC823 and 
HGC27 cells. Loss of expression was found in AGS cells 
and reduced expression was found in MGC803 and 
SGC7901cells (Figure 1A). The methylation status of 
the promoter region was determined by MSP. Complete 
methylation was found in AGS cells and partial meth- 
ylation was found in MGC803 and SGC7901 cells. No 
promoter region methylation was detected in NCI-N87, 
MKN45, BGC823 and HGC27 cells (Figure 1B). Lost/ 
reduced expression of HOXAII correlated with the 
promoter region methylation status very well. Resto- 
ration of HOXALI expression was induced by 5-AZA 
treatment in methylated cancer cells (Figure 1A). These 
results suggest that expression of HOXAII was con- 
trolled by promoter region methylation. To further vali- 
date the specificity of the MSP studies, and to determine 
the methylation density in the promoter region, sodium 
bisulfite sequencing was performed in AGS, BGC823 


and SGC7901 cells. As shown in Figure 1B, sequencing 
results matched the MSP results well, with near com- 
plete, partial or absent DNA methylation. 


HOXA11 was frequently methylated in human 
primary GCs & reduced expression correlates 
with promoter region hypermethylation 
HOXAI1 was found to be methylated in 81% (91/112) 
of primary gastric cancers, while no methylation was 
found in five cases of normal gastric mucosa (Figure 1C). 
Methylation of HOXALI was associated with male gen- 
der (p < 0.05), tumor size (p < 0.05), tumor differentia- 
tion (p < 0.05) and lymph node metastasis (p < 0.05), 
but no association was found with age and TNM stage 
(Table 1). To determine if there was an association of 
HOXALI expression and promoter region methylation 
in primary human GCs, the expression of HOXAI1 
was evaluated using immunohistochemistry in 45 cases 
using matched GCs and adjacent tissue. HOXA11 
staining was observed in both the nucleus and cyto- 
plasm. In 45 cases of GC samples, absent or reduced 
expression was found in 37 cases. In 45 cases of adja- 
cent tissue samples, absent or reduced expression was 
found in five cases. This suggests that the expression 
of HOXAII was significantly reduced in GC samples 
compared with normal tissues (p < 0.05, Figure 1D). 
However, of greatest importance was the correlation to 
promoter region DNA methylation. In the 37 cases of 
GC samples with absent or reduced levels of HOXAII 
expression, 36 cases were found to be methylated. In 
contrast, no methylation was found in the other eight 
cases of GC which express HOXA11. Thus repression 
of HOXAII was significantly associated with promoter 
region hypermethylation (p < 0.05, y? test, Figure 1D). 


HOXA11 suppresses cell proliferation in AGS 
cells 

To evaluate the effect of HOXALI expression in gastric 
carcinomas, cell viability and colony formation were 
used with cell lines with different baseline expression of 
HOXA11. The colony formation decreased after restora- 
tion of HOXAII expression in AGS cells 572 + 69 versus 
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Table 1 Clinicopathological features and HOXA11 methylation status in gastric cancer. 


Methylated 81.25% Unmethylated 


Clinical factor Total (n) 

(n = 91) 
Age (years): 
-<60 49 42 
->60 63 49 
Gender 
— Male 83 71 
— Female 29 20 
Alcohol abuse: 
— Negative 33 30 
— Positive 79 61 
Smoking: 
— Negative 37 31 
— Positive 75 60 
Tumor size (cm): 
-<5 37 25 
->25 75 66 
Differentiation: 
— Poor 61 56 
— High/moderate 51 35 
Tumor stage: 
-Il 16 12 
-IH 71 61 
-IV 25 18 
Lymph node metastasis: 
— Positive 88 76 
— Negative 24 15 


*p-values are obtained from a % test, significant difference: p < 0.05. 


HOXAT11 methylation status 


p-value (y? test) 
18.75% (n = 21) 


p = 0.2857 
7 
14 

p = 0.0490* 
12 
9 

p = 0.0905 
18 

p = 0.6294 
6 
15 

p = 0.0092* 
12 
9 

p = 0.0018 
5 
16 

p = 0.2431 
4 
10 
7 

p = 0.0183* 
12 
9 


116 + 30 (p < 0.05, Figure 2B). This difference suggests 
that colony formation was suppressed by HOXALI. This 
activity was further confirmed by shRNA knockdown 
in HOXALI expressing BGC823 cells. The colony for- 
mation increased from 77 + 8 to 281 + 21 (p < 0.05, 
Figure 2A) after shRNA knockdown of HOXAII. 

The cell viability was determined next using 
MTT. The OD value was 1.343 + 0.126 versus 
0.668 + 0.098 (p < 0.05) before and after restoration 
of HOXAII expression in AGS cells. The OD value 
was 2.536 + 0.342 versus 1.582 + 0.181 (p < 0.05) 
before and after knocking down HOXAII in BGC823 
cells (Figure 2B). These results suggest that HOXAI1 
suppresses cell proliferation in GC cells. 


HOXA11 induced G2/M arrest in GC cells 

To further analyze the effect of HOXAII on gastric 
carcinogenesis, cell phase distribution was evalu- 
ated by flow cytometry in AGS cells before and after 
re-expression of HOXAI1. As shown in Figure 2C, the 
distribution of cell phases in the AGS cell line with 
HOXAII unexpressed and re-expressed were: phase 


G0/1: 38.964 1.17% versus 34.78 + 3.57%; S phase: 
42.02 + 1.74% versus 27.13 + 5.80%; G2/M phase: 
19.68 + 2.34% versus 46.08 + 5.18%. The S phase was 
reduced significantly (p < 0.05) and the G2/M phase 
was increased significantly (p < 0.05). To further vali- 
date these findings, shRNA knockdown technique was 
used in HOXAII stably expressed BGC823 cells. The 
cell phase distribution, before and after knocking down 
HOXAII, was as follows: GO/Iphase: 57.22 + 0.07% 
versus 58.08 + 2.32%; S phase: 19.42 + 0.05% versus 
40.53 + 2.43%; G2/M phase: 22.98 + 0.59% versus 
1.38 + 0.14% (Figure 2C). S phase was increased signifi- 
cantly (p < 0.05) and the G2/M (p < 0.05) was reduced 
significantly after knocking down HOXAII. These 
results suggest that HOXAII induces G2/M phase 
arrest in GC. To further understand the mechanism, 
the expression of G2/M phase regulators, cyclinB1 and 
cdc2, was examined by western blot before and after 
re-expression of HOXAII in AGS cells. The expression 
of cyclin Bl and cdc2 was reduced after re-expression of 
HOXAII (Figure 2C), providing an explanation for the 
observed cell cycle changes. 
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Figure 2. The effect of HOXA11 on gastric cancer cell proliferation and apoptosis (see previous page). 

(A) Upper: Representative results of colony formation in HOXA17 unexpressed and re-expressed AGS cells. 
Lower: Representative results of colony formation in BGC823 cells before and after knocking down HOXA71. 
The experiment was repeated for three-times. *p < 0.05. (B) Upper: Growth curves in HOXA11 unexpressed and 
re-expressed AGS cells analyzed by MTT assay. Lower: Growth curves before and after knocking down HOXA11 
in BGC823 cells. Time Points: 6, 24, 48 and 72 h. Each experiment were repeated for three-times. *p < 0.05. 

(C) Left: Cell phase distribution evaluated by flow cytometry in AGS and BGC823 cells. Middle: Diagrams show the 
cell phase distribution of AGS and BGC823 cells. Right: The expression of cyclin B1 and cdc2 detected by western 
blot in HOXA11 unexpressed and expressed AGS cells. Actin: internal control. The experiment was repeated for 
three-times. *p < 0.05. (D) Left: Apoptosis results in HOXA17 unexpressed and expressed AGS cells. Right: The 
expression of caspase-3, survivin detected by western blot in HOXA11 unexpressed and expressed AGS cells. 


Actin: internal control. 


HOXAT1 induces apoptosis in human GC 

To further understand the role of HOXAII in GC, 
apoptosis was analyzed by flow cytometry. The rate of 
apoptotic cells was 7.41 + 0.151% versus 12.01 + 1.73% 
before and after re-expression of HOXAI/ in AGS cells 
(p < 0.05, Figure 2D). Apoptosis related proteins, cas- 
pase-3 and survivin were detected by western blot. The 
expression of caspase-3 was increased and the expression 
of survivin was reduced obviously after re-expression of 
HOXAII in AGS cells (Figure 2D). These results suggest 
that HOXAII induced apoptosis in human GC. 


HOXA11 inhibits GC cell migration & invasion. 
Since methylation of HOXAII is associated with 
lymph node metastasis in primary GC, we explored 
the effect of HOXAII on cell migration and inva- 
sion in GC cells. The number of migrating cells 
decreased from 250 + 12 to 98 + 8 (p < 0.05) after 
re-expression of HOXAII in AGS cells. The number 
of invasive cells was also decreased from 142 + 20 to 
63 + 10 (p < 0.05) after re-expression of HOXAII in 
AGS cells (Figure 3A). These effects of HOXAI were 
confirmed by knocking down HOXAII in BGC823 
cells. As shown in Figure 3A, the number of migrating 
cells increased from 242 + 34 to 827 + 24 (p < 0.05) 
cells and the number of invasive cells increased from 
170 + 15 to 356 + 21 (p < 0.05) after knocking down 
HOXAII in BGC823 cells (Figure 3B). 

The mechanism by which HOXAII altered cell 
migration and invasion was further explored by deter- 
mining the expression of MMP2 and MMP9, pro- 
teins associated with disruption of extracellular matrix 
which promote metastases. The expression of MMP2 
and MMP9 were both reduced after re-expression of 
HOXAII in AGS cells (Figure 3C). These results fur- 
ther support the role of HOXA/1 in suppression of cell 


migration and invasion in GC. 


Isolation and identification of cancer-related 
genes regulated by HOXA77 in GC 

To understand the mechanism by which HOXAII 
promotes carcinogenesis, gene expression micro- 
array was used. The expression of 223 genes was 


increased and 117 genes was reduced by more than 
twofold in HOXAII re-expressed AGS cells. Among 
these genes, 56 were upregulated more than four- 
fold and 21 were down regulated more than fourfold 
(Figure 4A, Supplementary Table 1, see online at www. 
futuremedicine.com/doi/suppl/10.2217/epi.14.92). 
These changes in gene expression were further vali- 
dated using RT-PCR. of greatest significance was the 
finding that WKDI was consistently up-regulated by 
re-expression of HOXAJI in AGS cells (Figure 4B). 
Combined with Pathway Association Analysis (http:// 
david.niaid.nih.gov), we identified NKDI, a key com- 
ponent of Wnt signaling for further study. Upregula- 
tion of NKDI was further determined by western blot. 
The effect of HOXAII up-regulating NKDI expres- 
sion was also confirmed in the knockdown of HOXA1I 
in BGC823 cells (Figure 4B). 


HOXAT1 inhibits Wnt/B-catenin signaling in GC 
cells 

To confirm that HOXAII involves Wnt signaling in 
GC development, TOPflash and FOPflash reporter 
systems were employed in this study. As shown in 
Figure 4C, the activity was reduced significantly by 
cotransfecting HOXAII with wild-type or mutant 
B-catenin with TOPflash reporter in AGS cells, but no 
changes were found by cotransfecting HOXALI with 
B-catenin and FOPflash reporter. 

The expression of Wnt signaling downstream genes 
was next examined using western blot. As shown in 
Figure 4D, the expression of c-Myc and cyclinD1 were 
both reduced and phosphorylated B-catenin (p-B- 
catenin) was increased after re-expression of HOXAI1 
in AGS cells, even though the total B-catenin level 
was not significantly changed (Figure 4D). The expres- 
sion of p-B-catenin was further evaluated in 45 cases 
of available matched primary human GC and adja- 
cent tissue samples by IHC. p-f-catenin staining was 
detected in both cytoplasm and nucleus. In 45 cases 
of adjacent tissue samples, 36 cases were p-P-catenin 
positive and nine cases were negative. In 45 cases of 
cancer tissue samples, six cases were p-B-catenin posi- 
tive and 39 cases were p-B-catenin negative (Figure 4E). 
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Figure 3. The effect of HOXA11 on gastric cancer cell migration and invasion. (A) Upper: Transwell assay shows migration results in 
HOXA11 unexpressed and re-expressed AGS cells. (200x) Middle: Transwell assay shows migration results before and after knocking 
down HOXA11 in BGC823 cells.(200x) Lower: Bar graphs represent the numbers of migrating AGS and BGC823cells. The experiment 


was repeated three-times. *p < 0.05. (B) Upper: Transwell assay shows invasion results in HOXA11 unexpressed and re-expressed 


AGS cells. (200x) Middle: Transwell assay shows invasion results before and after knocking down HOXA11 in BGC823 cells. (200x) 


Lower: Bar graphs represent the numbers of invasive AGS and BGC823cells. The experiment was repeated three-times. *p < 0.05. 


(C) The expression of MMP-2 and MMP-9 detected by western blot in HOXA11 unexpressed and re-expressed AGS cells. 


HOXA11 staining and p-P-catenin staining were asso- 
ciated (p < 0.05). These results suggest that HOXAII 
function as a Wnt signaling inhibitor in human GC. 


Discussion 

Homeobox sequences encode the 60-amino-acid 
homeodomain, with a helix-turn—helix DNA binding 
motif [25]. In mouse and man, there are approximately 
200 homeobox genes, with the 39 most homologous 
to the Drosophila homeotic genes located in four clus- 
ters (HOXA, HOXB, HOXC and HOXD) and termed 
Hox genes [26]. HOXAII is located in chromosome 
7p15.2, a region frequently deleted in human cancers 
[27-29]. Larger deletions of the HOXA gene cluster 
(HOXAI-HOXA13) result in hand-foot-genital syn- 
drome in combination with other malformations such 
as velopharyngeal insufficiency and a persistent patent 
ductus botalli [30]. A single HOXAII truncating muta- 
tion has been reported to cause a specific syndrome 


with skeletal defects and amegakaryocytic thrombo- 
cytopenia (31). HOXAII methylation was reported in 
different cancers, and methylation of HOXAII has 
been proposed as an early detection marker in breast 
cancer [32] and a poor prognostic marker in ovar- 
ian cancer [12]. In a Phase II clinical trial, decitabine 
demethylation of HOXAII restored sensitivity to 
platinum in patients with platinum-resistant ovarian 
cancer [33]. Bai et al. reported that HOXAII expression 
was silenced in human GC [34]. In this study, we found 
that HOXAII was frequently methylated in human 
GC and expression of HOXAII was regulated by pro- 
moter region methylation. Methylation of HOXAI1 
was associated with male gender, large tumor size, 
poor tumor differentiation and lymph node metasta- 
sis. HOXAII suppressed GC cell proliferation, inva- 
sion and migration and induced apoptosis. We found 
NKDI, a key component of Wnt signaling [35], was 
upregulated by HOXAII in GC cells, suggesting Wnt 
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Figure 4. HOXA11 is involved in Wnt signaling pathway in gastric cancer. (A) Gene-expression array results show differentially 
expressed genes in HOXA11 unexpressed and re-expressed AGS cells. Upper: The percentage of upregulated or down regulated 
genes were shown in the pie chart. Lower: Representative data of pathway-related genes according to pathway-related gene 
analysis [24]. (B) The expression of NKD1 was detected by semiquantitative RT-PCR and western blot in AGS and BGC823 cells. Actin: 
internal control. (C) Dual-Luciferase reporter assay. Luciferase activity was normalized to Renilla luciferase activity.*p < 0.05. (D) The 
expression of c-Myc, cyclinD1, p-B-catenin and B-catenin was detected by western blot in HOXA11 unexpressed and re-expressed AGS 
cells. Actin: internal control. (E) IHC results of p-B-catenin expression in adjacent normal tissue and gastric cancer tissue samples. 
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signaling as a pathway through which HOXAII medi- 
ates these effects. In mice, HOXAII has previously 
been shown to interact with Wnt signaling in regu- 
lating female reproductive tract development [36]. We 
demonstrate that HOXAII was indeed inhibits Wnt 
signaling as demonstrated with a Luciferase reporter 
assay in human GC cells. The inhibiting role of 
HOXAII in Wnt signaling was further confirmed by 
detecting the expression of Wnt signaling downstream 
genes in GC cells, as well as by evaluating HOXA11 
and p-B-catenin staining in human primary GC. 
Thus, HOXA11 suppresses GC growth and metastasis 
by inhibiting Wnt signaling. 

In primary tumors, quantitative methylation anal- 
ysis can be considered, however, all samples retain 
detectable unmethylated DNA which is derived from 
normal tissue contamination and unmethylated tumor 
cells in tumors with altered DNA methylation. Thus, 
distinguishing the degree of methylation present is 
often more a measure of tumor cellularity than any 
functional measure of silencing. For this reason, we 
have not distinguished a level of methylation, but 
compared the presence of methylation to expression 
to demonstrate functional consequences of this altera- 
tion. This may alternatively have been done by real- 
time MSP with the establishment of cutoffs. Some 
genes can become aberrantly methylated during cell 
culture. While, our purpose is to explore the conse- 
quences of DNA methylation on silencing and down- 
stream phenotypes in human GC. Promoter region 
methylation is found in 91 of the 112 primary GCs, 
suggesting that this is not solely a cell culture change. 

In addition to 5mC (5-methylcytosine), reduction 
in ShmC (5-hydroxymethylcytosine) was well reported 
in a number of cancers and cell lines corresponding to 
prostate, breast, colon, lung, brain, liver, kidney and 
melanoma, compared with the associated normal tissue 
[37-41]. But the level of ShmC is highly variable and un- 
reproducible [39,42]. It was undetectable in some human 
cancer tissues and cells [43]. ShmC of HOXA11 was 
not detected in this study. The mechanism of 5hmC in 
human GC in HOXAII remains need to be elucidated. 

The regulation of HOX gene expression has been 
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a model to study the establishment and maintenance 
of heritable transcriptional states during development 
[44]. Polycomb (PcG) and trithorax (TrxG) proteins 
maintain appropriate HOX expression. The PRC2 
complex trimethylates H3K27(me3) [45] and bind to 
the same histone mark that is essential for inheritable 
long-term repression of target genes [46]. Recruitment 
of PRC1 members to H3K27me3 further facilitates 
gene silencing. While, TrxG/MLL protein complexes 
catalyze the trimethylation of H3K4, which is gener- 
ally associated with active transcription [47]. The regu- 
lation of HOXA11 expression by histone modification 
remain unclear. To better understand the mechanism 
of HOXAII expression regulating, it is necessary to 
analyze the promoter region histone modification 
status in GC. 


Conclusion 

HOXA11 is frequently methylated in human GC and 
the expression of HOXA1I was regulated by promoter 
region methylation. Restoration of HOXAII expres- 
sion suppressed cell proliferation, migration and inva- 
sion, as well as induced apoptosis and G2/M arrest. 
HOXAII suppresses GC growth by inhibiting Wnt 
signaling. HOXAII methylation is a potential GC 
detection and prediction marker. HOXAII is a Wnt 
signaling inhibitor that could become a therapeutic 
target in GC. 


Future perspective 

HOXAI1 methylation is a potential gastric can- 
cer detection and prediction marker, prognosis. 
As HOXAII is a Wnt signaling inhibitor, it is possible 
become a gastric cancer therapeutic target. 


Financial & competing interests disclosure 

This work was supported by grants from the National Ba- 
sic Research Program (973 Program No. 2012CB934002, 
2010CB912802), National Key Scientific instrument Special 
Programme of China (Grant No. 2011YQ03013405), National 
High-tech R&D Program (863 Program No. $S2012AA020314, 
$S2012AA020821, SS2012AA020303) and National Sci- 
ence Foundation of China (Grant No. 81121004, 81071953, 


HOXA11 inhibiting Wnt signaling by upregulating NKD1 expression in human gastric cancer 
HOXA11 is frequently methylated in human gastric cancer (GC) and the expression of HOXA77 was regulated 


by promoter region methylation. 


Promoter region methylation of HOXA11 was associated with male gender, tumor size, tumor differentiation 


and lymph node metastasis. 


HOXA71 suppressed GC cell proliferation, invasion and migration, as well as induced apoptosis and G2/M 


phase arrest. 


NKD1, a key component of Wnt signaling, was up-regulated by HOXA717. 


HOXAT1 suppresses GC by inhibiting Wnt signaling. 


fsg 


future science group 


www.futuremedicine.com 


211 


Research Article Cui, Gao, Linghu et al. 


81161120432). The authors have no other relevant affiliations 
or financial involvement with any organization or entity with a 
financial interest in or financial conflict with the subject mat- 
ter or materials discussed in the manuscript apart from those 
disclosed. 

No writing assistance was utilized in the production of this 
manuscript. 


References 


1 Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin 
DM. Estimates of worldwide burden of cancer in 2008: 
GLOBOCAN 2008. Int. J. Cancer 127(12), 2893-2917 
(2010). 


2. Jemal A, Bray F, Center MM, Ferlay J, Ward E, Forman D. 
Global cancer statistics. CA Cancer J. Clin. 61(2), 69-90 
(2011). 


3 Kelley JR, Duggan JM. Gastric cancer epidemiology and risk 
factors. J. Clin. Epidemiol. 56(1), 1-9 (2003). 


4 Milne AN, Carneiro F, O’morain C, Offerhaus GJ. Nature 
meets nurture: molecular genetics of gastric cancer. Hum. 


Genet. 126(5), 615—628 (2009). 


5 Panani AD. Cytogenetic and molecular aspects of gastric 
cancer: clinical implications. Cancer Lett. 266(2), 99-115 
(2008). 


6 Ali Z, Deng Y, Tang Y, Zheng S, Ma N, He N. 
Epigenetic deregulations in gastric cancer. J. Nanosci. 
Nanotechnol. 13(1), 40—51 (2013). 


7 Park JH, Park J, Choi JK et al. Identification of DNA 
methylation changes associated with human gastric cancer. 
BMC Med. Genomics 4, 82 (2011). 


8 Duboule D, Dolle P. The structural and functional 
organization of the murine HOX gene family resembles that 
of Drosophila homeotic genes. EMBO J. 8(5), 1497-1505 
(1989). 


9 Mcginnis W, Krumlauf R. Homeobox genes and axial 
patterning. Cell 68(2), 283-302 (1992). 


10 Gross S, Krause Y, Wuelling M, Vortkamp A. Hoxall and 
Hoxd11 regulate chondrocyte differentiation upstream of 
Runx2 and Shox2 in mice. PLoS ONE 7(8), e43553 (2012). 


11. Cheng W, Liu J, Yoshida H, Rosen D, Naora H. Lineage 
infidelity of epithelial ovarian cancers is controlled by HOX 
genes that specify regional identity in the reproductive tract. 
Nat. Med. 11(5), 531-537 (2005). 

12. Fiegl H, Windbichler G, Mueller-Holzner E et al. HOXAII1 
DNA methylation — a novel prognostic biomarker in ovarian 
cancer. Int. J. Cancer 123(3), 725-729 (2008). 

13 Whitcomb BP, Mutch DG, Herzog TJ, Rader JS, Gibb RK, 
Goodfellow PJ. Frequent HOXAII and THBS2 promoter 
methylation, and a methylator phenotype in endometrial 
adenocarcinoma. Clin. Cancer Res. 9(6), 2277-2287 (2003). 

14 Hwang JA, Lee BB, Kim Y et al. HOXAII hypermethylation 
is associated with progression of non-small cell lung cancer. 
Oncotarget 4(12), 2317-2325 (2013). 

15 Apostolidou S, Hadwin R, Burnell M et al. DNA 


methylation analysis in liquid-based cytology for cervical 
cancer screening. Int. J. Cancer 125(12), 2995-3002 (2009). 


Ethical conduct 

The authors state that they have obtained appropriate institu- 
tional review board approval or have followed the principles 
outlined in the Declaration of Helsinki for all human or animal 
experimental investigations. In addition, for investigations in- 
volving human subjects, informed consent has been obtained 
from the participants involved. 


16 Yan W, Wu K, Herman JG et al. Epigenetic silencing of 
DACHI induces the invasion and metastasis of gastric cancer 
by activating TGF-beta signalling. J. Cell. Mol. Med. 18(12), 
2499-2511 (2014). 


17. Yan W, Wu K, Herman JG et al. Epigenetic regulation of 
DACHI, a novel Wnt signaling component in colorectal 
cancer. Epigenetics 8(12), 1373-1383 (2013). 


18 Herman JG, Graff JR, Myohanen S, Nelkin BD, Baylin SB. 
Methylation-specific PCR: a novel PCR assay for methylation 
status of CpG islands. Proc. Natl Acad. Sci. USA 93(18), 
9821-9826 (1996). 


19 Jia Y, Yang Y, Liu S, Herman JG, Lu F, Guo M. SOX17 
antagonizes WNT/beta-catenin signaling pathway in 
hepatocellular carcinoma. Epigenetics 5(8), 743-749 (2010). 


20 Derks S, Lentjes MH, Hellebrekers DM, De Bruine AP, 
Herman JG, Van Engeland M. Methylation-specific PCR 
unraveled. Cell Oncol. 26(5—6), 291-299 (2004). 


21 Cregger M, Berger AJ, Rimm DL. Immunohistochemistry 
and quantitative analysis of protein expression. Arch. Pathol. 
Lab. Med. 130(7), 1026-1030 (2006). 


22 Koo CL, Kok LF, Lee MY et al. Scoring mechanisms 
of p16INK4a immunohistochemistry based on either 
independent nucleic stain or mixed cytoplasmic with nucleic 
expression can significantly signal to distinguish between 
endocervical and endometrial adenocarcinomas in a tissue 
microarray study. J. Transl. Med. 7, 25 (2009). 


23 Jia Y, Yang Y, Brock MV, Zhan Q, Herman JG, Guo 
M. Epigenetic regulation of DACT2, a key component 
of the Wnt signalling pathway in human lung cancer. 
J. Pathol. 230(2), 194-204 (2013). 


24 Huang DW, Sherman BT, Lempicki RA. Systematic 
and integrative analysis of large gene lists using DAVID 
Bioinformatics Resources. Nat. Protoc. 4(1), 44—57 (2009). 


25 Qian YQ, Billeter M, Otting G, Muller M, Gehring 
WJ, Wuthrich K. The structure of the Antennapedia 
homeodomain determined by NMR spectroscopy in 
solution: comparison with prokaryotic repressors. Cell 59(3), 
573-580 (1989). 


26 Zhao Y, Potter SS. Functional comparison of the Hoxa 4, 
Hoxa 10, and Hoxa 11 homeoboxes. Dev. Biol. 244(1), 21-36 
(2002). 


27 Sander S, Bullinger L, Leupolt E et al. Genomic aberrations 
in mantle cell lymphoma detected by interphase fluorescence 
in situ hybridization. Incidence and clinicopathological 
correlations. Haematologica 93, 680—687 (2008). 

28 Nucci MR, Harburger D, Koontz J, Dal Cin P, Sklar J. 
Molecular analysis of the /AZFI—JJAZI gene fusion by RT- 


PCR and fluorescence in situ hybridization in endometrial 
stromal neoplasms. Am. J. Surg. Pathol. 31(1), 65-70 (2007). 


212 


Epigenomics (2015) 7(2) 


fsg 


future science group 


The mechanism of HOXA717 in gastric carcinogenesis 


Research Article 


29 Prokunina-Olsson L, Fu Y-P, Tang W er al. Refining the 38 Jin SG, Jiang Y, Qiu R er al. 5-Hydroxymethylcytosine 
prostate cancer genetic association within the JAZFI gene is strongly depleted in human cancers but its levels do 
on chromosome 7p15.2. Cancer Epidemiol. Biomarkers not correlate with JDH/ mutations. Cancer Res. 71(24), 

Prev. 19(5), 1349-1355 (2010). 7360-7365 (2011). 

30 Devriendt K, Jaeken J, Matthijs G et al. Haploinsufficiency 39 Laird A, Thomson JP, Harrison DJ, Meehan RR. 
of the HOXA gene cluster, in a patient with hand-foot- 5-hydroxymethylcytosine profiling as an indicator of cellular 
genital syndrome, velopharyngeal insufficiency, and state. Epigenomics 5(6), 655—669 (2013). 
persistent patent ductus Botalli. Am. J. Hum. Genet. 65(1), 40 Nestor CE, Ottaviano R, Reddington J et al. Tissue type is 
249-251 (1999). a major modifier of the 5-hydroxymethylcytosine content of 

31. Fryssira H, Makrythanasis P, Kattamis A et al. Severe human genes. Genome Res. 22(3), 467—477 (2012). 
developmental delay in a patient with 7p21.1-p14.3 41 Yang H, Liu Y, Bai F et al. Tumor development is associated 
microdeletion spanning the TWIST gene and the HOXA with decrease of TET gene expression and 5-methylcytosine 
gene cluster. Mol. Syndromol. 2(1), 45—49 (2011). hydroxylation. Oncogene 32(5), 663-669 (2013). 

32 Tommasi S, Karm DL, Wu X, Yen Y, Pfeifer GP. 42 Ficz G, Gribben JG. Loss of 5-hydroxymethylcytosine in 
Methylation of homeobox genes is a frequent and early cancer: cause or consequence? Genomics 104(5), 352-357 
epigenetic event in breast cancer. Breast Cancer Res. 11(1), (2014). 

R14 (2009). i 
43 Dahl C, Gronbaek K, Guldberg P. Advances in DNA 

33 Matei D, Fang F, Shen C er al. Epigenetic resensitization to methylation: 5-hydroxymethylcytosine revisited. Clin. Chim. 
platinum in ovarian cancer. Cancer Res. 72(9), 2197-2205 Acta 412(11-12), 831-836 (2011). 

(2012); 44 Soshnikova N, Duboule D. Epigenetic regulation 

34 Bai Y, Fang N, Gu T et al. HOXAII gene is hypermethylation of vertebrate Hox genes: a dynamic equilibrium. 
and aberrant expression in gastric cancer. Cancer Cell Epigenetics 4(8), 537-540 (2009). 

Int. 14(1), 79 (2014). i i 
45 Cao R, Wang L, Wang H et al. Role of histone H3 lysine 27 

35 Van Raay TJ, Coffey RJ, Solnica-Krezel L. Zebrafish Naked methylation in Polycomb-group silencing. Science 298 (5595), 
and Naked2 antagonize both canonical and non-canonical 1039-1043 (2002). 

War signaling. Dev: Biol 3092), 151-1682007). 46 Francis NJ, Follmer NE, Simon MD, Aghia G, Butler JD. 

36 Mericskay M, Kitajewski J, Sassoon D. Wnt5a is required Polycomb proteins remain bound to chromatin and DNA 
for proper epithelial-mesenchymal interactions in the uterus. during DNA replication in vitro. Cell 137(1), 110-122 
Development 131(9), 2061-2072 (2004). (2009). 

37 Haffner MC, Chaux A, Meeker AK er al. Global 47 Briggs SD, Bryk M, Strahl BD et al. Histone H3 lysine 4 
5-hydroxymethylcytosine content is significantly reduced methylation is mediated by Setl and required for cell growth 
in tissue stem/progenitor cell compartments and in human and rDNA silencing in Saccharomyces cerevisiae. Genes 
cancers. Oncotarget 2(8), 627-637 (2011). Dev. 15(24), 3286-3295 (2001). 

fsg www.futuremedicine.com 213 


future science group 


